Capstan friction model for DNA ejection from bacteriophages 
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Bacteriophages infect cells by attaching to the outer membrane and injecting their DNA into 
the cell. The phage DNA is then transcribed by the cell's transcription machinery. A number 
of physical mechanisms by which DNA can be translocated from the phage capsid into the cell 
have been identified. A fast ejection driven by the elastic and electrostatic potential energy of the 
compacted DNA within the viral capsid appears to be used by most phages, at least to initiate 
infection. In recent in vitro experiments, the speed of DNA translocation from a A phage capsid 
has been measured as a function of ejected length over the entire duration of the event. Here a 
mechanical model is proposed that is able to explain the observed dependence of exit velocity on 
ejected length, and that is also consistent with the accepted picture of the geometric arrangement 
of DNA within the viral capsid. 



The physics of semi-flexible polymers under confine- 
ment is a subject of great interest with wide applications 
in biology and soft matter. Viruses are an interesting ex- 
ample of such an application. Bacteriophages are viruses 
that infect bacteria. During infection, the virus capsid 
itself remains outside the bacterium, only the DNA is in- 
jected into the cell pQ. Several physical mechanisms by 
which bacteriophages inject DNA into their hosts have 
been identified. The specific mechanism varies by phage 
type and stage of the infection process [2H1] • A fast ejec- 
tion on the timescale of seconds can be achieved simply 
by the release of elastic and electrostatic energy of the 
coiled up DNA confined within the viral capsid. It has 
been shown in in vitro studies that A phage can eject 
its entire 48.5 kbp genome by this mechanism in a few 
seconds [3D]. In other phage species, this "syringe" mech- 
anism may initiate the infection, which is then completed 
through other pathways [5J. 

A schematic diagram of the DNA ejection process is 
shown in Fig. fl| If F(N) denotes the free energy of the 
compacted DNA, then —F'(N)/b is the force driving the 
ejection, where N is the number of base pairs of DNA 
confined within the capsid and b w 0.31 nm is the dis- 
tance between base pairs. The free energy of the part of 
the DNA that is already outside the viral capsid is ne- 
glected as it is much smaller than the contribution from 
the confined part. Since inertia is completely negligible 
at such small scales, this driving force must be balanced 
by a frictional force, vm(N), where v is the DNA translo- 
cation velocity and m is a frictional coefficient (mobility) 
which may in general depend on N. The linear depen- 
dence on the velocity is a consequence of the fact that 
at low Reynolds numbers, the fluid equations reduce to 
the linear Stokes equation. Thus, if the functions F(N) 
and m{N) were known, the translocation velocity could 
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FIG. 1: Sketch showing the arrangement of the viral DNA 
in the capsid. The right panel shows a sectional view when 
the point of observation is on the central axis. The DNA 
enters the capsid through the tail opening (central circle on 
the right panel) and is forced to curve on encountering the 
wall. It then wraps around the central axis in helical coils of 
decreasing radius. 

be determined from 

-F'(N)/b = vm(N). (1) 

The function F(N) can be calculated analytically [71 H] 
on the presumption that the DNA is packed in the capsid 
as helical coils wound in successive layers of decreasing 
radii in an "inverse spool" arrangement (Fig. [T]). It has 
been shown that such an ordered packing does indeed 
minimize the free energy of the system. In Brownian 
dynamics simulations [9] , as the chain is gradually intro- 
duced into the capsid, it is observed to arrange itself in 
a donut shape. As N becomes larger, it winds into a 
spool from the outside in. When the internal diameter 
of the spool shrinks sufficiently, the DNA stops wrapping 
in helical coils and simply forms loose turns parallel to 
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the spool axis. The most direct evidence for the helical 
conformation of DNA within the capsid has come from 
imaging of DNA inside T7 [TO] and T4 [TT] phages by 
cryo-electron microscopy and of P22 and A phages by X- 
ray diffraction [TJ]. The images reveal multiple layers of 
helical coils of DNA wound around the central axis. 

The analytical formula for the free energy F(N) has 
been subjected to direct experimental tests by a rather 
ingenious method. Real cells exert significant osmotic 
pressures that are believed to resist DNA insertion by 
the phage. It has been possible to actually "stall" the 
DNA translocation in in vitro experiments by raising the 
osmotic pressure in the solution outside the capsid. In 
such cases, the translocation stops with an amount N = 
N n of DNA still remaining within the capsid, where Nq 
is given by 



F'(N )/b = A Pos 



(2) 



Here Ap os is the difference in osmotic pressures between 
the outside and inside the capsid and a is the DNA ra- 
dius. The quantities Ap os and N can be directly mea- 
sured for a range of values: from Ap os — to a value that 
is high enough that none of the DNA is ejected. Thus, 
F(N) can be determined experimentally and compared 
to analytical models based on the elastic and electrostatic 
energy of a helically coiled semi-flexible rod. This has 
been done [T31 HI] and the theoretical models for F(N) 
were found to agree very well with experimental data. 
Further evidence for the model comes from studying an 
event that occurs much later in the infection process: the 
packaging of viral DNA into the capsid prior to lysis of 
the host cell and release of the phages into the environ- 
ment. In an in vitro experiment, Smith et al [15] attached 
an optical bead to a DNA strand being actively packaged 
into the 029 bacteriophage capsid by molecular motors. 
By applying a variable force to the optical bead and ob- 
serving the changes in the packing rate, they concluded 
that the elastic force against which the molecular motors 
must work, increases rapidly with the amount of DNA in- 
side the capsid. This is again consistent with the rapidly 
decreasing nature of the function F(N). 

Unfortunately, no such theoretical result for the func- 
tion m(N) in equation (I]) is available. However, if v 
is measured experimentally, then using equation (JlJ the 
function m(N) can be obtained. This was done in a re- 
cent experiment by Grayson et al. |16j where A phages 
were immobilized on a surface and induced to eject their 
DNA by exposing them to cell surface proteins to which 
the phage would normally attach and trigger the ejec- 
tion event in the native state. The ejected DNA was 
visualized by a fluorescent dye and the ejection process 
recorded by video microscopy. The data from the exper- 
iment (Fig. 5 in their paper) is reproduced in Fig. [2] fro m 
which certain conclusions may be drawn about the ori- 
gin of the frictional resistance encountered by the DNA 
during the ejection process. 

Frictional resistance to the motion of DNA can arise 
either from friction within the phage tail or from move- 
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length of DNA in capsid, N [kbp] 

FIG. 2: Dependence of mobility on DNA length remaining in 
capsid. Measured values of m(N) from Fig. 5 in the paper by 
Grayson et al. [16] shown with the fitted exponential, equa- 
tion djjb. The data are for AcI60 in NaCl buffer (solid circle), 
Ab221 in NaCl buffer (solid square), AcI60 in MgS0 4 buffer 
(open circle) and Ab221 in MgS04 buffer (open square). Error 
bars are not shown. 



ment of the DNA outside or within the capsid. In the 
former case, m(N) would be independent of N. Fig. [2] 
shows that m(N) becomes approximately independent 
of N only towards the very end of the translocation pro- 
cess after 80 to 90 percent of the DNA has already been 
ejected. This is consistent with the following simple esti- 
mate of the frictional drag that is obtained if the DNA is 
taken as a cylinder of radius a ~ 1 nm, concentric with a 
slightly larger cylinder (the phage tail) of radius a* ~ 1.5 



/drag M 



-2naL 



(3) 



where L ~ 150 nm is the length of the tail and /i « 10~ 3 
Pa-s is the viscosity of water. If we take v ~ 60 kbp/s, 
the highest velocity indicated by the data, then /drag ~ 
0.04 pN may be regarded as the maximum frictional drag 
arising from the tail. This is much less than the 10 — 40 
pN driving force [H] ejecting the DNA from the capsid 
and therefore cannot be the principal source of resistance. 
Similarly, it is easy to show that viscous friction from 
the medium external to the capsid cannot be important 
either. The actual drag would depend on the shape of 
the ejected DNA coil, but, an estimate, that should be 
regarded as an upper bound, is obtained by regarding 
the entire N max =48.5 kbp DNA to be a cylinder that 
translates in water parallel to its axis at a speed of v ~ 60 
kbp/s. In this case, the drag force is [17] 



/dra 



ln(&AW»+ln2-l/2 



0.4 pN, (4) 



which is again negligible in comparison to the driving 
force. However, the driving force does decay rapidly with 
N, and for N < 5 kbp the driving force is less than a pN. 
Thus, the hydrodynamic drag could balanced the driving 
force towards the very end of the translocation process, 
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and indeed, the data shows that m(N) stops decreasing 
once N becomes smaller than about 10 kbp (Fig. [2]). 

Fig. [2] also shows that as long as N > 10 — 20 kbp, the 
experimental data can be fit very well by a function of 
the form 



/drag/w = m(N) = Aexp(fciV). 



(5) 



The fitted straight line corresponds to A — 0.012 pN- 
s/kbp and k = 0.105 kbp -1 . It is significant that m(N) 
is hardly affected when the buffer contains the divalent 
Mg 2+ ion which reduces the electrostatic repulsion be- 
tween the DNA strands thereby reducing the driving 
force by almost an order of magnitude. This is consis- 
tent with our expectations, since frictional resistance is 
not expected to be affected greatly by interstrand repul- 
sion. 

The success of the empirical fit, equation ([5]), makes 
one wonder if there is any known physical mechanism 
that could lead to such an exponential friction law. We 
show here that if one assumes that the resistance to mo- 
tion arises due to friction between the sliding DNA strand 
and its neighboring strands as well as with the capsid 
wall, then such an exponential law can be derived. For 
simplicity, we consider the shape of the viral capsid to 
be a cylinder, though bacteriophage capsids are typi- 
cally polyhedral in shape and are sometimes modeled as 
spheres. We start from the equations of equilibrium of 
an elastic beam [15] 



dF 

ds 
dM 

ds 



— L 
F x t 



(6) 
(7) 



where s is the arc length along the beam, F is the elastic 
force across a cross-section, M is the internal bending 
moment and L is the external force per unit length along 
the beam. Since the pitch of the helix is small [7], we will 
neglect its torsion. The external force is, by Amonton's 
law [19], L = —Nh — fNt, where t and n are the unit 
tangent and normal to the curve representing the beam 
centerline, t points in the direction of sliding and / is a 
friction coefficient. The force across the cross section may 
be resolved into a tensile (T) and shear («S) component: 
F = Tt + Sh. The bending moment is related to the local 
curvature (k) as M = Elnh, where b = t x n, E is the 
Young's modulus and / is the area moment of inertia of 
the cross-section. Substitution in equations Q and ([T]) 
and use of the Frenet-Serret formulas J20j describing how 
the unit vectors (t,n, b) change along a curve in three 
dimensions results in 
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(8) 
(9) 
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For a helix, k is a constant, thus S = 0, and therefore, 
dT 



ds 



which may be integrated to yield 

Tx = T o exp(f0) 



(11) 



(12) 



where T\ is the tension in the DNA at the capsid exit 
and To is the tension at the terminal point. By terminal 
point we mean the point on the DNA within the capsid 
marking the transition from a helix to the loosely coiled 
trailing end in the central core. If we express the viscous 
resistance on this trailing end as vm , then T — vtuq. 
Since the pitch of the helix is small, <f> s» 2nbN/ R, so that 



m(N) = m exp(27r& fN/R), 



(13) 



where R ~ 25 nm, is the capsid radius. It should be 



noted that in equation (12), Ti/Tq is independent of the 



radius of the turns but depends only on the total turn 
angle, <fi. Thus, if there are multiple layers of coils, it may 
still be applied as long as the friction coefficient between 
all contacting surfaces are the same. Equation ( 13 ) is 
equivalent to equation ([5]) ; comparing the coefficients wc 
determine / ~ 1.5. 



It is interesting to note that equation ( 12 ) is the cap- 



stan equation (also known as Eytelwein's formula or the 
Euler-Eytelwein formula) familiar from mechanics [21) . 
In the classical "capstan problem" a flexible line is wound 
around a cylinder (for example, a bollard or capstan used 
in ship mooring) to hold a load at one end of the line by 
the application of a holding force. At the point when 
the line is just about to slip, the relation between the 
load (Ti) and holding force (To) is described by equa- 
tion (12). The exponential dependence of T\/Tq on 



explains how a very large load can be supported by a rel- 
atively modest force if the rope is wound even a few turns 
around the capstan. The name "capstan friction model" 
in the title of this paper is a reference to this correspon- 
dence. It should be noted, however, in our case there is 
no central cylinder. The frictional resistance arises from 
contact between neighboring DNA strands and possibly 
with the capsid wall as the tightly coiled DNA unwinds. 
In fact, the unfurling of a surveyor's steel measuring tape 
is perhaps a better analogy. 

The frictional coefficient between a pair of surfaces de- 
pends on the nature of the contact between them [T91I22], 
When a lubricating fluid film is present, the coefficient 
is small but increases linearly with sliding speed in a 
manner that may be well understood from the equa- 
tions of fluid flow. In the presence of a large normal 
load (which in the present problem would arise from the 
strong bending rigidity of the DNA) the fluid film thins 
and one enters the regime of boundary lubrication where 
direct molecular level contact between surfaces could oc- 
cur. This is characterized by a sharp rise in the frictional 
coefficient. The value / ~ 1.5 corresponds to this latter 
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regime of boundary lubrication. For example, the fric- 
tional coefficient between dry highly polished metal sur- 
faces at moderately large loads is close to this value [TH] • 
The derivation of equation ( 13 ) presumes Amonton's law 



of friction which is based on the idea that contact be- 
tween a pair of surfaces arise from interlocking asperities 
and is thus independent of the apparent area of contact 
but increases in proportion to the normal load |19) . Such 
a picture is unlikely to hold for nanoscale phenomena 
such as the sliding of DNA strands. Nevertheless, Amon- 
ton's law is frequently used to describe frictional effects in 
nanoscale phenomena even though the underlying phys- 
ical models are different from the rubbing of asperities 
envisaged to rationalize Amonton's law in the classical 
treatment of friction [23TE5] . 

It may be shown by simple estimates that viscous fric- 
tion with the fluid in the core is unimportant. Since any 
shear induced at the outer boundary of the fluid filled 
core will equilibrate on a time scale R 2 Jv ~ 1 ns [y 
is the kinematic viscosity of water), we conclude, that, 
the fluid core must be in rigid rotation during most of 
the ejection process of duration r ~ 1 s. The inertia 
of this fluid mass also plays no role. Indeed, the mo- 
ment of inertia of the cylinder is / ~ pR 5 and the rate 
of change of angular speed is w ~ bN max /Rr 2 where 
N m ax is the total length of DNA. Substituting typical 
values JTS] we find that ICj ~ 10~ 32 N-m which is many 



orders of magnitude smaller than the applied torque, 
M RF'(N)/b ~ 1CT 18 N-m. Thus, neither the vis- 
cosity of the water within the capsid nor its inertia plays 
a significant role in determining the ejection velocity. 

In this paper we interpret phenomena involving the 
translocation of single molecules of DNA across nanome- 
ter size pores from the principles of classical continuum 
mechanics. In the light of the small scale nature of 
the system, one might question the validity of such an 
approach. However, even though the relevant length 
scales are approaching the limits of applicability of the 
continuum description of matter, classical continuum 
mechanics has been applied successfully to other prob- 
lems involving translocation of single molecules across 
nanopores Thus, its use in this context, for 

estimating forces on DNA, appears entirely reasonable. 
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